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Abstract. Until recently, the physically relevant singularities occurring in FRW cosmologies 
had traditionally been thought to be limited to the big bang, and possibly a big crunch. 
However, over the last few years, the zoo of cosmological singularities considered in the 
literature has become considerably more extensive, with big rips and sudden singularities 
added to the mix, as well as renewed interest in non-singular cosmological events such as 
bounces and turnarounds. In this talk, we present an extensive catalogue of such cosmological 
milestones, both at the kinematical and dynamical level. First, using generalized power 
series, purely kinematical definitions of these cosmological events are provided in terms of 
the behaviour of the scale factor a(t). The notion of a scale- factor singularity is defined, 
and its relation to curvature singularities (polynomial and differential) is explored. Second, 
dynamical information is extracted by using the Friedmann equations (without assuming even 
the existence of any equation of state) to place constraints on whether or not the classical 
energy conditions are satisfied at the cosmological milestones. Since the classification is 
extremely general, and modulo certain technical assumptions complete, the corresponding 
results are to a high degree model- independent. 



1. Introduction 

The physically relevant singularities in cosmology have traditionally been thought to be limited 
to the big bang or big crunch ,l a . However, over the last few years, more singularities have been 
considered such as big rips, sudden singularities, and there has even been a renewed interest 
in non-singular events such as bounces and turnarounds (21 El IU El El • 

From observational cosmology we extract the cosmological principle, which states that the 
universe is homogeneous and isotropic. Hence, we consider a cosmological spacetime of the 
FRW form: 

ds 2 = -di 2 + a(t) 2 | - + r 2 [dO 2 + sin 2 9 d0 2 ] j , (1) 

where a(t) is the scale factor of the universe and k = ±1,0 is the space curvature. 

Assuming the applicability of the Einstein equations of general relativity, we can then deduce 



the Priedmann equations. In geometric units where 8ttGn = 1 and c = 1, we have 



a' cr 

. . a a 2 
pi) = -2 2--, 3 

p(t)+3p(t) = -6~, (4) 
a 

and the related conservation equation 

p(t) a 3 + 3 \p(t) + p(i)] a 2 a = 0, (5) 

where any two of equations (0 El EI) imply ® . 

In this talk we shall derive generic definitions of all the singularities considered above (which 
we shall refer to as cosmological milestones) and we shall pin down necessary and sufficient 
conditions for their occurence both at the level of kinematics and dynamics. 

• For the kinematical approach, we develop precise definitions of these cosmological 
milestones in terms of the scale factor a(t), and with the notion of scale factor singularity 
in hand we check whether these milestones are true curvature singularities. 

• For the dynamical approach, we use the Friedmann equations to analyze the classical 
energy conditions in the vicinity of these milestones in a model-independent manner. 

2. Kinematical analysis 

2.1. Cosmological milestones definitions 

From a mathematical point of view, solutions of a problem can often be expanded in Taylor 
series or Laurent series around their singular points. Hence, the general idea is to expand the 
scale factor a(t) in generalized power series, similar to a Puisieux series, in the vicinity of the 
cosmological milestones. 

Generic cosmological milestone: Suppose we have some unspecified generic cosmological 
milestone, that is defined in terms of the behaviour of the scale factor a(t), and which occurs 
at some finite time t@ . We will assume that in the vicinity of the milestone the scale factor has 
a (possibly one-sided) generalized power series expansion of the form 

a{t) = co\t - tol* +d\t- i©!" 1 + c 2 |i - t & \^ +c 3 \t- top + . . . (6) 

where the indicial exponents rji are generically real (and are often non-integer) and without 
loss of generality are ordered in such a way that they satisfy 

m< vi <m<m-- ■ ( 7 ) 

Finally we can also without loss of generality set 

c > 0. (8) 
There are no a priori constraints on the signs of the other q, though by definition q ^ 0. 



The first term of the right hand side of equation © is the dominant term, and is therefore 
responsible for the convergence or divergence of the scale factor at the time £q. The indices rji 
are used to classify the cosmological milestones and the absolute values are used to distinguish 
a past event from a future event. This power series expansion of the scale factor is sufficient 
to represent almost all the physical models that we are aware of in the literature. Table ^ 
represents this cosmological milestone classification depending on the value of the scale factor. 
For Big bangs or Big crunches the scale factor is zero and it corresponds to 770 > in terms of 



Table 1. Classification of cosmological milestones 



Cosmological 


Scale factor 


Indices 


milestones 


value 


Vi 


Big Bang/ 


o(t ) = 


Vo > 


Big Crunch 






Sudden 


a(t®) = c 


Vo = 


Singularity 


a (")(t ) = 00 


7/1 non-integer 


Extremality 


a(t Q ) = c 


Vo = 


events 




Vi € Z+ 


Big rip 


a(i Q ) = 00 


Vo < 



the power series expansion indices. The value 770 = leads to the scale factor being finite and 
describes two types of milestones, sudden singularities and extremality events — which are not 
singularities in any sense (bounce, turnaround, inflexion). Specifically, sudden singularities are 
of order n where the n th derivative of the scale factor is infinite: 

a (n) (t t Q ) ~ c vi (Vi ~ l)(m -2)...( m -n + l)\t- t Q \^~ n 00, (9) 

and therefore 771 has to be a non- integer [IJ[S|. In contrast for extremality events, the power 
series expansion of the scale factor a(t) is a simple Taylor series with r\i E Z + . Furthermore, 
an extremality event is of order n where: 

• a( 2n )(i©) > for a bounce (local minimum) [7||B], 

• a( 2n )(i©) < for a turnaround (local maximum), 

so that for a bounce or turnaround we write: 

a(t) = a(t G ) + -I-a^)(i ) [t - t & ] 2n + ...; (10) 

• a^ 2n+1 \tQ) 7^ for an inflexion event (neither local minimum nor maximum), 
so that: 

a(t) = a(t Q ) + ^-l—^a( 2n+1 Ht Q ) [t - t G ] 2n+1 + .... (11) 

Finally, a big rip occurs when the scale factor a(t) is infinite at some finite time t© |H El, it 
corresponds to 770 < in terms of the power series expansion. Even though we distinguish a 
"future rip" from a "past rip" , the literature to date has only considered future rips as a past 
rip would be a very unusual beginning to the universe. 



Summary: We have written a generic definition of the scale factor a(t) based on generalized 
power series (generalized Puisieux series) for all the physically relevant cosmological milestones 
found in the literature to date (big bang, big crunch, sudden singularity, extremality events and 
big rip). We can now use the parameters to explore the kinematical and dynamical properties 
of the cosmological milestones, see whether they are true curvature singularities and whether 
the energy conditions hold in the vicinity of the time of the event t©. Note that for most 
calculations it is sufficient to use the first three terms of the power series expansion. 



2.2. Spacetime curvature 

To determine whether a cosmological milestone is a true spacetime curvature singularity or 
not, in the sense of classical general relativity, we need to have a closer look at the Riemann 
tensor at the time of the event. However, due to the symmetries of the FRW geometry we are 
assuming, the Weyl tensor vanishes and hence it is sufficient to look at the Ricci tensor. Now 
if we take into account the spherical symmetry and translational symmetry, we only need to 
look at the only two non-zero orthonormal independent components of the Ricci tensor to test 
for curvature singularities: 

Rff and Rff = R§§ = R^- (12) 

On another hand, one could consider instead Rff and the Ricci scalar R, or Rff and the Einstein 
tt-component Gff, or alternatively the Ricci scalar R and the combination R a i,R ab . 



Polynomial curvature singularity: 

Proposition 1. To characterize all polynomial curvature singularities in a FRW spacetime 
geometry, and hence decide whether the cosmological milestones we have defined are true 
curvature singularities, it is sufficient to test Rff and Gff in orthonormal components for 
finiteness: 

% = - 3 | <%-s(£ + *). ,13) 

Note that Rff is independent of the space curvature k, Gff is independent of the second 
derivative of the scale factor d, and finally Rff and Gff are linearly independent. 

From the definition of a generic cosmological milestone (assuming t > t & for simplicity) we 
have: 



o(t) = c (t-t Q )^ + Cl (t-t )'n + ... (14) 
a(t) = c r ]0 (t-t Q r- 1 + c m (t-t Q T 1 - 1 + ... (15) 
a(t) = c Qm { m -l)(t-t Q )^- 2 + c m ( m -l)(t-t & )^- 2 

+c2V2(v2-m-t G r- 2 + ... (i6) 

The method we use to study the finiteness of Rff or Gff is generic. We write Rff in the vicinity of 
a cosmological milestone using the power series expansion of the scale factor and its derivatives, 
keeping only the most dominant terms: 

• if rjo ^ and rji ^ 1 then 

R tt = -3 3 ^ < f } °, S , hence lim R u = sign(r/ [l - %]) oo (17) 

a It — if) r t^tQ 



• if 770 = and r/i ^ 1 then 
771(171 - l)ci 



% ~ "3 



co 



(t-t ) ,?1 - 2 , hence lim % 



= < 





-6ci/c 
sign(ci)oo 



771 > 2; 
m = 2; 
m G (1,2); 



-sign(ci)oo 771 G (0, 1); 

(18) 



• if both r/o = and 771 = 1 then 

i? t - t - ~ -3 m ^ 2 - 1)02 (t - i©)" 2 " 2 , hence lim % 







co 



-6c 2 /c 



% >2; 
r/2 = 2; 



-sign(c 2 )oo 772 G (1,2); 

(19) 



• if 770 = 1, then (since 771 > 1) 
771(771 - l)ci 



co 



(i-i ) m ~ 3 , hence lim % 

t— >*0 







771 > 3; 
-18ci/c 771 = 3; 

-sign(ci)oo 771 G (1,3); 



(20) 

Now we can see that the number of cases in which remains finite is rather limited. Explicitly, 
i?££ is finite provided that: 

• Vo = 0, m > 2 ; 

• 770 = 0, 771 = 1, 772 > 2; 

• 770 = 1 and 771 > 3. 

We apply the same generic argument to analyze G^ in the vicinity of a cosmological milestone. 
We have: 



• if 770 7^ then 

• while if 770 = then 



G it ~ 3 



Vo 



+ 



G u ~ 3 



[{t-t & f cl(t-t & )^ 



(21) 



(22) 



After a few technical calculations we find that the necessary and sufficient conditions for 
to remain finite are: 

• Vo = 0, 771 > 1; 

• 770 = 1, Cq + k = 0, and 771 > 3. 

(Note that Cq + k = implies k = —1, and cq = 1.) 

To decide the characterization of a cosmological milestone as a polynomial curvature singularity, 
we need to combine the two previous results. 

Proposition 2. Assuming we can write the scale factor as a generalized power series in a the 
vicinity of a cosmological milestone, the necessary and sufficient conditions for both G^ and 
to remain finite, so that this cosmological milestone is not a polynomial curvature singularity, 
are: 



• Vo = 0, m > 2; 

• m = 0, 771 = 1, 772 > 2; 

• 770 = 1, k = —1, Co = 1, and 7/1 > 3. 

Note that the first two cases concern a sub-class of sudden singularities (7/0 = 0) and the 
last case corresponds to a type of big bang or big crunch (7/0 = 1). This result seems quite 
surprising, however, one can ask whether those exceptional cases are excluded if we generalize 
the idea of polynomial curvature singularity to that of derivative curvature singularity. 

Derivative curvature singularity: Even if the Riemann tensor does not blow up in the vicinity 
of a cosmological milestone, does this property still hold for some derivative of the Riemann 
tensor? We define a n^-order derivative curvature singularity by the condition that the n 
derivative of the curvature tensor blows up. Because of the symetries of the FRW universe that 
we mentioned earlier the only derivatives of interest are time derivatives, and hence we only 
need to consider: 

a (n+2) 

—3 h (lower-order derivatives); (23) 

a a( n+1 ) 

3 k h (lower-order derivatives). (24) 

a 1 

Proposition 3. To avoid a n ttl -order derivative curvature singularity in the vicinity of a 

(j) 

cosmological milestone, we must force all - — with j ^ n + 2 to remain finite. Subsequently, 
the indices rji in the generalized power series expansion of the scale factor a(t) must be positive 
integers forcing the series to be a Taylor series. 

Conclusion: Are cosmological milestones true curvature singularities? Assuming we can write 
a(t) as a generalized power series in the vicinity of some cosmological milestone, we have 
demonstrated that, with the exception of two very specific cases, almost all these events are 
true curvature singularities in the sense that they are both polynomial and derivative curvature 
singularities. The two situations in which a cosmological milestone is not a true curvature 
singularity are if: 

• Vo = 0, rji £ Z + ; corresponding to an extremality event (bounce, turnaround, or inflexion 
event) rather than a bang, crunch, rip, or sudden singularity; 

• 7/o = l, k = — 1, Co = 1, t]i € Z + , and rji > 3; corresponding to a FRW geometry that 
smoothly asymptotes near the cosmological milestone to the Riemann-flat Milne universe 1 
which is not a really popular model in cosmology. 

3. Dynamical analysis 

We now wish to explore the consequences of the previous kinematical results at a dynamical 
level using the Friedmann equations ((21 El EI an d © ■ How "strange" does physics get in the 
vicinity of some cosmological milestone? To answer this question we ask what happens to the 
energy conditions of general relativity at the cosmological milestones. Even though one can 
question if they are truly fundamental 9 , they help quantifying how "strange" physics gets at 
the cosmological events by informing us as to which type of physical conditions are satisfied or 
not. 



d n R 



it 



d n t 
d n t 



1 The case a(t) — t, k = 



— 1 is the Milne universe and is a disguised portion of Minkowski space. 



The energy conditions of general relativity are the null, weak, strong, and dominant energy 
conditions. Assuming a FRW universe and making use of the Friedmann equations (0 and 
EJ), they specialize to the following conditions 003113: 

• The [NEC] is satisfied if p+p > that is, in view of the Friedmann equations, if k > a a— a 2 . 

• The [WEC] is satisfied if the NEC is and if in addition p > 0, that is k > —a 2 . 

• The [SEC] is satisfied if the NEC is and if in addition p + 3p > 0, that is a < 0. 

• The [DEC] is satisfied if p ± p > 0, which reduces to the NEC and the condition 
k > -(a a + 2d 2 )/2. 

Note that the NEC is the most interesting condition as it is the weakest of all and it leads 
to the strongest theorems (e.g. Hawking area theorem for black hole horizon). The DEC is 
usually satisfied for "normal matter" and can be interpreted as saying that the speed of energy 
flow of matter is always less than the speed of light. 

We now consider the generalized power series expansion for the scale factor a(t) defined 
in (|14|) and the consecutive derivatives ()15j) and ()16|) . Following the same pattern detailed in 
section 12.21 for Rii, we obtain a complete classification of all the energy conditions as being 
satisfied or not depending on the different values of the space curvature k, the indices rji, and 
the coefficients Cj (see ^U] for full details). In this article we will not present the classification 
in full detail, but we will mention some of the main results for the different energy conditions. 

NEC: The classification obtained for the NEC is in agreement with most of the results 
scattered throughout the literature. For example, it is obvious that big rips always violate 
the NEC and therefore all the other energy conditions. Note that the NEC holds for big bangs 
and big crunches independently of the space curvature k within the range r/o G (0, 1). However 
for a class of sufficiently "violent" bangs (r/Q > 1) and hyperbolic spatial curvature (k = — 1) 
the NEC is violated. The analysis becomes more tedious for the classes of sudden singularities; 
many but not all of them violate the NEC. 

WEC: Remember that the condition k > —a 2 , in addition to the NEC, must hold in order 
for the WEC to be satisfied. Consequently, the WEC is equivalent to the NEC concerning 
bangs, crunches and rips, and the only changes arise for the class of cosmological milestone 
corresponding to r/o = (sudden singularities or extremality events). The WEC also holds for 
big bangs and big crunches independently of the space curvature k within the range r]Q G (0, 1). 
Again note that many but not all of the sudden singularities violate the NEC. 

SEC: Not only do the big rips violate the NEC, WEC, and therefore the SEC, but all "violent" 
bangs or crunches violate the SEC (independently of k). However, the SEC is satisfied for big 
bangs and big crunches independently of the space curvature k within the range 770 G (0,1). 
Also, note that turnarounds satisfy the SEC while bounces violate the SEC. Regarding inflexion 
events, the SEC is violated either just before or just after (in agreement with [3IH]). 

DEC: There are several classes of cosmological milestones for which the DEC holds all the 
way to the singularity. In particular, the DEC is satisfied for big bangs and big crunches 
independently of the space curvature k, within the more restrained range t/q G (1/3, 1). More 
surprisingly, there are a few classes of sudden singularities satisfying the DEC. This result 
seems to be in contradiction with Lake's analysis jS], however the disagreement resides in the 
definition of just "how sudden" is a sudden singularity. Only some sufficiently "gentle" sudden 



singularities, those for which the first divergence occur in higher derivatives a^ n \t — > i Q ) for 
n > 3, satisfy the DEC. Indeed, for some of these events the DEC can be satisfied all the way 
to the singularity. This result is in agreement with Barrow and Tsagas jSj, though now we 
have a full characterization of those situations. 

4. Discussion 

Assuming a FRW universe, and applicability of general relativity, we have explored three issues 
concerning singularities: 

• We have developed a generic classification of the various physically relevant "cosmological 
milestones" occurring in the literature in terms of a generalized power series expansion of 
the scale factor a(t). 

• With the notion of generalized power series, we have defined (on a purely kinematical level), 
the idea of polynomial and derivative curvature singularities. It is therefore possible to 
classify all "cosmological milestones" as to whether they are singular or not, with the class 
of non-singular milestones being strictly limited (to extremality events or asymptotically 
Milne bangs/crunches). 

• On dynamical grounds, we have checked the validity or otherwise of the classical energy 
conditions in the vicinity of the cosmological milestones. Depending on one's attitude 
towards the energy conditions 9 , one could use this catalogue as a guide towards deciding 
on potentially interesting scenarios to investigate. 

Finally, if in the vicinity of any cosmological milestone, the input scale factor a(t) is a 
generalized power series, then all physical observables (e.g. H, q, the Riemann tensor, etc.) 
will likewise be a generalized power series. By checking the related indicial exponents that can 
be calculated from the indicial exponents of the scale factor, one can determine whether or not 
the particular physical observable then diverges at the cosmological milestone. 
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